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ABSTRACT 

Doty,  William,  K.,  An  Introduction  to  the  Application  of  Linear 
Programming  to  Industrial  Engineering,  May,  1954..  There  are  50  pages, 
5  figures,  and  10  references  in  the  bibliography. 

The  purpose  of  this  paper  is  to  explain  what  linear  programming 
is,  how  the  linear  programming  problem  is  interpreted  mathematically, 
how  to  solve  a  simple  problem,  what  the  assumptions,  limitations,  and 
possibilities  of  this  technique  are,  and  how  to  set  up  for  solution  a 
few  examples  of  feasible  problems.  The  emphasis  is  on  Industrial 
Engineering  applications  and  requires  no  knowledge  of  mathematics  be- 
yond college  algebra. 

The  mathematical  discussion  illustrates  the  necessary  elementary 
vector  algebra  including  addition  of  vectors,  unit  vectors,  and  linearly 
independent  sets.  Linear  transformations  are  illustrated  by  transfor- 
mations from  a  3-dimensional  to  a  2-dimensional  space.  Convex  sets 
and  convex  polyhedral  cones  are  illustrated.  A  simple  numerical  example 
with  three  unknowns  in  two  equations, with  a  linear  functional  to  be 
maximized,  is  solved  step  by  step  and  the  geometrical  significance  of 
each  step  is  explained. 

The  four  assumptions  necessary  for  the  application  of  linear 
programming  are  discussed  and  certain  feasible  problem  areas  in 
Industrial  Engineering  are  listed.  It  is  shown  why  the  technique  can- 
not be  used  for  scheduling,  sequence  of  production  selection,  elimina- 
tion of  expediters,  etc. 

Three  mathematical  models  are  presented:  (1)  machine  assignment 
to  maximize  total  productivity,  (2)  production  line  loading  for  a  single 


product  to  meet  sales  requirements  over  successive  time  periods  and  to 
minimize  the  combined  costs  of  regular  and  overtime  production  and 
storage,  and  (3)  a  numerical  example  to  show  the  complexity  to  which 
even  problems  of  small  size  are  subject.  This  model  deals  with  three 
products,  two  production  operations,  three  time  periods,  and  the  max- 
imization of  profit. 

The  conclusion  points  out  that  the  application  of  linear  program- 
ming to  Industrial  Engineering  is,  at  present,  limited,  but  predicts 
that,  as  study  in  both  theory  and  practical  aspects  continues,  more 
and  more  Industrial  Engineering  application  will  be  made  in  the  future. 


AN  INTRODUCTION  TO  THE  APPLICATION  OF 
LINEAR  PROGRAMMING  TO  INDUSTRIAL  ENGINEERING 

INTRODUCTION 
Considerable  interest  has  developed  in  the  past  few  years  in 
the  mathematical  technique  known  as  linear  programming.  This  interest, 
at  first,  centered  itself  among  mathematicians  and  economists  because 
the  concepts  were  similar  to  or  were  advancements  in  technique  under 
study  in  those  two  sciences.  Interest  in  other  fields  was  limited 
until  Dr.  George  B.  Dantzig  evolved  the  "simplex  method"  of  computa- 
tion. This  was  a  turning  point.  The  formerly  theoretical  and  limited 

aspects  of  linear  programming  were  now  found  to  have  direct  application 

2         3  / 

to  practical  problems  facing  the  armed  forces  ,  industry  ,  agriculture^ 

and  transportation  .  Evidence  of  widening  interest  is  illustrated  by 

the  appearance  of  articles  on  linear  programming  in  such  practical 

minded  publications  as  Factory  Management  &  Maintenance  and  Business 


1.  Dantzig,  G.3.,  "Maximization  of  a  Linear  Function  of  Variables 
Subject  to  Linear  Inequalities,"  Activity  Analysis  of  Production 
and  Allocation,  Cowles  Commission  Monograph  13;  New  York,  John  Wiley 
&  Sons,  1951.  PP  339-3A7. 

2.  Wood,  M.K.,  and  Geisler,  M.A.,  "Development  of  Dynamic  Models  for 
Program  Planning,"  (ibid.),  pp.  189-215. 

3.  Charnes,  A.,  Cooper,  W.W.,  and  Mellon,  B.,  "Blending  Aviation 
Gasolines,"  Econometrica.  vol.  2Ct   p.  2. 

4.  Hildreth,  C,  and  Reiter,  S.,  "On  the  Choice  of  a  Crop  Rotation 
Plan,"  Activity  Analysis  of  Production  and  Allocation,  Cowles 
Commission  Monograph  13;  New  York,  John  Wiley  &  Sons,  1951,  pp.  177- 
188 

5.  Dantzig,  G.B.,  "Application  of  the  Simplex  Method  to  a  Transporta- 
tion Froblem,"  (ibid.),  pp.  359-373. 

6.  Melden,  M.E.,  "Operations  Research,"  Factory  Management  &  Mainten- 
ance, October  1953,  vol.  Ill,  No.  10,  pp.  113-120. 


Week  .  In  addition,  the  success  of  Operations  Research  and  the  otten- 

/  '» 
dent  publicity  has  helped  to  focus  attention  to  linear  programing 

because  this  mathematical  technique  is  a  useful  tool  for  Operations 

Research  teams . 


7.  "Top  Management  by  Mathematics,"  Business  Week,  May  30,  1953, 
No.  1239,  pp.  86-92. 


PURPOSE 

The  majority  of  the  literature  to  date  has  been  slanted  to  the 
field  of  economics  and  mathematics  and  hence  requires  a  certain  degree 
of  mathematical  maturity  for  its  understanding.  A  few  articles 
appearing  in  periodicals  have  intimated  the  possibilities  of  applica- 
tion of  linear  programming  to  Industrial  Engineering  problems.  Between 
the  two  approaches  lies  a  wide  gap. 

The  purpose  of  this  thesis  then  is  to  attempt  to  bridge  this  gap 
by  explaining  what  linear  programming  is,  how  the  linear  programming 
problem  is  interpreted  mathematically,  how  to  solve  a  simple  problem, 
what  the  assumptions,  limitations,  and  possibilities,  of  this  technique 
are,  and  how  to  set  up  for  solution  a  few  examples  of  feasible  problems. 
The  present  treatment  will  endeavor  to  cover  the  material  without 
requiring  any  mathematics  beyond  college  algebra. 


PRELIMINARY  DISCUSSION 

Linear  programming  is  a  mathematical  technique  which  may  be 
utilized  in  the  optimal  planning  of  activities  which  are  interdepen- 
dent. Programming  in  the  title  indicates  the  area  of  interest,  the 
preparation  of  plans,  which,  when  executed,  will  yield  the  best 
results.  Programming  is  apparent  in  such  areas,  for  instance,  as 
mixtures  of  chemicals,  loading  of  machine  tools,  assignment  of  person- 
nel, and  shipping  allocation. 

By  linear  is  meant  the  restriction  that  the  variables  or 
unknowns  must  occur  to  the  first  power.  No  squares,  cubes,  etc.,  are 
permissable,  nor  may  one  variable  be  multiplied  by  another.  Best  or 
optimal  results  may  refer,  for  example,  to  maximizing  profits,  or 
production,  or  to  minimizing  cost  of  production  or  storage. 

The  linear  programming  technique  is  carried  out  as  follows: 

1.  A  mathematical  model  (or  set  of  equations)  is 
formulated  from  the  word  problem,  this  set  of 

equations  to  follow  a  certain  form. 

2.  The  mathematical  model  is  solved , usually  by  the 
simplex  method. 

In  regard  to  the  formulation  of  the  mathematical  model  it  need 
not  apply  that  the  equations  in  the  model  express  rigorously  every 
facet,  every  fine  point  which  could  conceivably  affect  the  problem 
at  hand.  Assumptions  and  approximations  may  be  necessary.   "Standard 
time",  a  statistical  approximation,  for  example,  may  be  employed  in 
linear  programming.  Simplification  of  the  problem  may  be  advantageous, 
if  not  necessary.  A  certain  phase  of  a  production  activity,  for  exam- 
ple, may  be  ignored  if  its  effect  on  the  total  is  small.  When  formu- 


lating  the  model,  therefore,  assumptions,  approximations  and  simpli- 
fications may  be  made  provided  the  model  retains  enough  predictive 
value  to  be  useful. 

In  the  construction  of  the  mathematical  model  there  are  two 
methods  of  approach.  In  the  majority  of  the  literature  which,  as 
has  been  stated,  has  evolved  from  an  economics  background  where  the 

interest  lay  in  covering  an  entire  economy,  large  formalized  models 

g 
of  economic  activities  have  evolved  .  This  activity  analysis  ap- 
proach, as  it  is  called,  has  been  extended  to  the  theory  of  the  firm 
so  that  many  concrete  and  every  day  problems  of  industry  may  be 
investigated". 

The  second  method  of  approach  to  the  construction  of  the  model, 
which  is  the  one  advanced  by  Charnes,  Cooper  and  Henderson  ,  is  "to 
approach  each  problem  squarely,  forming  the  model  and  calculational 
procedures  naturally  from  the  demands  of  the  material."  They  feel  that 
the  activity  analysis  approach  is  awkward  and  artificial  and  "carries 
with  it  the  danger  of  forcing  the  problem  or  overlooking  important 
aspects  which  do  not  fit  easily  into  prepared  rubrics." 

The  construction  of  a  mathematical  model,  especially  for  a 
complex  problem,  however,  is  an  art  and  requires  a  keen  mind  coupled 
with  a  good  mathematical  and  technological  background .  Examples  of 
models  will  be  presented  later. 


8.  Leontief,  W.U.,  The  Structure  oT  the  American  Economy  1919-1939* 
second  ed.;  New  York,  Oxford  University  Fress,  1951. 

9.  Dorfman,  R.,  Application  of  Linear  Irogranmlng  to  the  Theory  of  the 
Firm;  Berkeley;  University  of  California  Fress,  1951. 

10.  Charnes,  A.,  Cooper,  W.W.,  and  Henderson,  A.,  An  Introduction  to 
Linear  Programming;  New  York,  John  Viley  &  Sons,  1953,  Preface. 


MATHEMATICS  OF  LINEAR  PROGRAMMING 
Before  illustrating  the  application  of  linear  programming  to  problems 
it  is  necessary  that  the  reader  understand  certain  basic  material  of 
vector  geometry  and  matrix  algebra.  Only  the  minimum  of  information 
necessary  to  the  solution  of  linear  programming  problems  will  be 
presented.  Mathematical  proofs  and  amplifying  material  can  be  found 
in  the  material  listed  in  the  bibliography,  especially  in  Part  II  of 

An  Introduction  to  Linear  Programming  by  Charnes,  Cooper  and  Hender- 

11 
son. 

Points  and  Vectors 

The  representation  of  a  point  in  a  plane  by  (x_,  y.)  where 
x.,  is  the  x  coordinate  (distance  from  the  origin  or  zero  point  along 
the  x  axis) ,  and  y,  is  the  coordinate  in  the  y  direction,  is  well 
known.  Just  as  (xj,  y^)  f    (x?>  y  )  locate  points  in  a  plane, 
(x  ,  y  ,  z.)  and  (x2,  y0,  z2)  locate  points  in  the  3-dimensional  space, 
(x,,  y-,  z-.),  in  addition  to  representing  a  point,  may  represent  a 
vector  which  can  be  visualized  as  an  arrow  whose  tail  is  at  the  origin 
and  whose  head  is  at  the  point  (x  ,  y. ,  z,). 

Similarly,  though  less  easily  visualized,  (x_  ,  x?,  x~,....x  ) 

represents  a  point  or  vector  in  an  n-dimensional  space,  where  x,,  x2, 

through  x  are  real  numbers  called  the  coordinates  in  the  dimensions 
n 

1,  2,  through  n.  Point  or  vector  x  may  be  represented  either  as 
1»  x2»,#,,»  xn) , which  is  called  a  row  vector, 

xl 
x2 

or  as      *       ,  which  is  called  a  column  vector. 


xn 


11.  Ibid. 


The  subscript  under  x  designates  the  coordinate  meant.  A 
superscript  above  identifies  to  which  point  one  refers.  Thus  x  and 
x  means  point  1  and  point  2  respectively  while  xi  means  the  second 
coordinate  of  point  3. 

Two  or  more  vectors  in  the  same  n-dimensional  space  may  be 
added  by  adding  corresponding  coordinates  to  form  a  new  vector.  Thus, 

if  x  =  (x£,  x2,  x3)  and  x  -  (xi,  x2,  x?)  then  x3,  the  sum  of  x1  and 

x  is  presented  as  x3  »  (x3,  x3,  x^)  .  (*J  *  xj,  x*  4  x2,  x*  +  x2)  with 


3    1,23    1.2       jl 
xi  s  x^  4  x_ ;  xl  s  x„  4  x  ;  x  -  x„  4-  x 

x    -1    1'  2    2    2 


3    1  ,  2 
3  =  X3  *  X3. 


This  addition  of  vectors  may  be  illustrated  in  2  dimensions  as 
follows:  Given:  x1  =  (2,1);  x2  =  (1,  3).  Find  x3,  the  sum  of  x1  and 
x2. 

X3  s  x1  4.  x2 

=  (xj  4  x2,  x*  4.  x2) 

=  (241,  143) 
=  (3,  4). 

The  point  x  '  may  be  found  also  by  laying  off  frcw  point  x^  '  a 

vector  parallel  to  and  the  same  length  as  vector  xv   ,  as  shown  in 
Figure  1. 


Fig.  1  Addition  of  Tvo  Vectors 


A  point  or  vector  may  be  multiplied  by  a  real  number.  Thus  if 
x  s  (2,  1),  then  2x  =   [(2)  (2)  ,   (2)(1)J   -  U,  2).  Such  a  multi- 
plication of  a  vector  by  a  real  number  gives  a  new  vector  whose  direc- 
tion is  exactly  the  same  as  the  original  vector  but  whose  length  is 
equal  to  the  length  of  the  original  vector  multiplied  by  the  number. 

A  unit  vector  is  a  vector  of  value  unity  (one)  in  one  axis  only, 
and  having  value  zero  in  every  other  axis.  They  are  designated  by 

e   t   e   »  e   ,  e   ,  etc.,  where  the  superscript  designates  which 

coordinate  has  value  one.  Thus  in  2-dimensional  space:  e^   =  (1,  0) 

(2) 

and  ev  '  =  (0,  1)  as  shown  in  Figure  2. 


2"    *? 


A  e 


(2) 


e 


(1) 


Fig.   2     Unit  Vectors  in  2-Dimensional  Space 

Similarly  e^   '    in  an  n-dimensional  space  is  represented  as: 
e^1)   =  (0,  0,..,  0,   1,  0,..,  0)   where  the  1  occurs  in  the  ith  dimension. 

Note   that  every  vector  in  a   2-dimensional  space   (the  plane)    can 
be  formed  by  a  linear  combination  of  the  two  unit  vectors  e^   '   and 
e'2'.     For  example:     the  vector  rS   '   -  (3,   2)    can  be  obtained  from 
(3)   e(l>   +  (2)   e<2)   ■  3   (1,  0)    ♦  2  (0,   1)   =  (3,   2). 

Similarly,   every  vector  in  an  s-dimensional  space  can  be  formed 
by  a  suitable  linear  combination  of  the  s  unit  vectors. 

These  s  unit  vectors  in  an  s-dimensional  space  together  form 

12 
what  is  known  as  a  linearly  independent  set.         This  means  that  for  a 

linear  combination  of  the  s  unit  vectors,  a^  e'    '    +  a2  e'^    +...-♦•  a     e^    , 

to  equal  the  zero  vector   (0,  C, ,  0),   the   coefficients  of  the 

e's  must  each  equal  zero.     In  2  dimensions  for  example: 

al  e(1)    *  a2  e(2)   =  (0>  0) 
r  ax   (l,  C)   ♦  a2  (C,   l)   =  (C,  0) 

=  (*1,  0)    4-  (0,  a2)   =  (0,  0) 


12.  Ibid.   p.  4-8. 
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=  (ai>  a2)  =  (0,  0)  only  if  ax  =  0  and  a?  =  0.  e^  and  e^2) 
are  thus  a  linearly  independent  set. 

No  vector  in  a  linearly  independent  set  can  be  formed  from  a  linear 
combination  of  the  other  vectors  in  the  set.  For  example,  no  linear 
combination  of  e^  '   and  e^  '  can  yield  e   . 

This  set  of  unit  vectors  is  not  the  only  linearly  independent 
set.  For  example,  in  2-dimensional  space,  (3,  1)  and  (2,  1)  are 
linearly  independent,  for  if  a(?,  1)  4-  b(?,  1)  m   (0,  0) ,  then 

1.  3a  4-  2b  =  0  and 

2.  a  +  b  s  0 

From  (2),  b  =  -a,  and  substituting  in  (1),  3a  -  2a  ;0  or 
a  =  0  and  thus  b  =  0. 

Note,,  however,  that  in  2-dimensional  space,  (1,  3)  and  (2,  6) 
are  not  linearly  independent  but  are  linearly  dependent.  Each  vector 
is  a  real  number  multiple  of  the  other:  2(1,  3)  r  (2,  6). 

An  additional  consequence  of  linearly  independent  sets  is  that 

if  a  point  P  r  a^1  +  a^2  4-  *  a^  (  where  x1,  x2, ,  xn  is 

a  linearly  independent  set),  then  the  coefficients  a  ,  a„, ,  a 

are  the  only  coefficients  of  the  x  's  in  terms  of  which  one  can  express 
P  as  a  sum  of  multiples  of  the  x  's,  i.e.,  as  a  linear  combination. 

The  maximum  number  of  points   (vectors)    in  a  linearly  independent 
set  in  2-dimensional  space   is   2,   in  3-dimensional  space  is  3,  and  in 
m-dimensional  space   is  m.     In  an  m-dimensional  space,   any  linearly 
independent  set  of  m  points  is  called  a  basis.     There  can  be  any  number 
of  bases  but  each  basis  of  an  m-dimensional  space  must  consist  of 
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exactly  m  linearly  independent  points. 

Linear  Transformations 
An  important  concept  in  linear  programming  is  that  of  linear 
transformations.  A  transformation  is  a  mapping  or  transfer  of  points 
from  an  n-dimensional  space  into  a  space  of  possibly  different  dimen- 
sions, each  point  assuming  its  new  position  in  accordance  with  the 
directions  specified  by  the  particular  transformation.  For  example, 
if  the  point  x  =  (1,  3)  were  transformed  by  T  to  the  point  (?,  6)  the 
transformation  T  of  x  equals  (2)  x  indicated  by  T  (x)  a  2x. 

A  linear  transformation  is  merely  a  restricted  type  of  trans- 

13 
formation  and  must  conform  to  the  following  two  requirements: 

1.  T(ex)  must  equal  aT(x) .  For  example,  if  x  -  (1,  3),  a  =  2, 
and  T  specifies  a  multiplication  by  3 

T(ax)  =  (3)(2)(1,  3)  Z   (6,  13) 

aT(x)  -  (2)(3)(1,  3)  =  (6,  18)  T  thus  conforms  to 
the  first  requirement. 

2.  T(x  *■   w)  must  eoual  T(x)  +  T(w)  . 

If  T,  again,  specifies  a  multiplication  by  3, 

x  s  (1,  3)  and  w  =  (2,  /,) ,  then 

T(x  t  w)  m   30  +2,  3*4)-  (S  21) 

T(x)  =  3(1,  3)  =  (3,  9)  and 

T(w)  =  3(2,  4)  =  (6,  1?)  adding, 

T(x)  *  T(w)  .  (3  +  6,  9  +  12)  ■  (9,  21) 

T  thus  conforms  to  the  second  requirement. 


13.  Ibid.  p.  U. 
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Because  the  set  of  ni  unit  vectors  in  an  m-dimensional  space  is 
a  linearly  independent  set,  or  basis,  every  point  in  the  m-dimensional 
space  can  be  expressed  as  a  linear  combination  of  the  m  unit  vectors. 
Thus,  if  x  :  (3,  0,  2)  ,  x  can  be  expressed  as 

x  =  (3)  e^  +  (0)  e(?)  +  (?)  e^)  or 

3 
x  =  \  ci  e'1'  where  c,  a  3;  c~  s  0;  c-  r  2 

i=  1 

It  follows   therefore   that  if  it  is  known  what  the  linear  trans- 
formation does  to  the  unit  vectors,   the  linear  transformation  of  any 
point  is  known. 

If  e'1'    b  (0,...,    1,  0,..,  0)   where  the  1  occurs   in  the  ith 

dimension,   and  i  =1,    ?, ,   n,   are  unit  vectors  in  an  n-dimensional 

space  N,   and  g^'    =  (0,..,1,  0,....,   0),  with  j   =  1,    ?,...,  m,  are  unit 

vectors  in  the  m-dimensional  space  M,   then  the   linear  transformation  T 

m 

of  e^  takes  the  form   T   e^   =   V"  ai .  g^J'^  .  Thus:  ^ 

T     ew     a.ugu)  ♦.128<:?)/..*.tagW,7(« 

T    e<«    =aklgW+ak2g<^t..  *vgW=yW 

The  coefficients  of  the  g^'*s  above  when  written  by  themselves 
can  be  written  in  an  array  as  follows: 


H.  Ibid.  p.  17. 
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all 

ai2 

• 

• 

Blm 

*21 

a?? 

• 

• 

aan 

akl 

\2 

• 

• 

km 

anl 

a 

• 

• 

a 
ran 

Such  an  array  is  customarily  called  a  matrix  of  n  rows  and  m  columns. 
The  matrix  is  said  to  be  associated   (^^^  )   with  the  linear  transform- 
at ion  T.     The  entries  in  the  k       row  are  the  coordinates  of  the  point 
into  which  e^    '    is  carried  by  the  transformation.     Thus,  e^   '   which  is 
the  point   (1,  C,  0,...,  0)    in  the  n-dimensional  space  N  is  transformed 
into  the  point  y^    '   whose  coordinates  are   (flh-i,   a-.^****  arm)    in  the 
m-dimensional  space  M. 
Example : 

N  is  a  3-dimensional  space  with 

e(1)   =  (1,  0,  0);  e(2)   =  (0,   1,  0);  e(3)   =  (0,  0,  1) 

M  is  a  2-dimensional  space  with 

g(1)    =  (1,  0);      g(?)   =  (0,   1),  and: 
2         1 

T      ^^         l  2 

2         3 

Then  T       e(1)       =  2  g(1)     +     1  g(r)     =     y1     »     (2,1) 

T       e(2)       =  1  g(1)     ♦     2  g(?)     =     y*     =     (1,   ?) 
T       e(3)       -  2  g^     ♦     3  g{2)     =    y^     -     (2,  3) 

Figure  3  illustrates  these  transformations  from  space  N  into 
space  M. 


u 


T(x) 


UA) 


y(3)=Te(3L(2,3) 


=  Te(2)=(l,?) 


y(1)=Te(1)=(2,l) 


-I h- 

U  5 


Fig.   3       Geometric  Representation  of  Linear  Trans formations 

Since  it  is  known  what  the  linear  trans  formation  T  does   to  the 
unit  vectors  in  H,  what  T  does  to  each  point  x  is  also  known.     If 
x  «  (xlf   x2,   x3)      «     (1,   0,   1), 
T   (x)      =     x     T  e^      4-     x2  T  e^2^      +     x^  T  e^ 

=     (1).(2,   1)      ♦     (0)    (1,    2)      +■     (1)    (2,   3) 

=     (2,1)      +     (2,3)      =     [(2  4-  2),    (1  +  3)]     =  U,  ti 

Points  x  and  T  (x)  are  illustrated  in  Figure  3 . 


Convex  Sets 
A  convex  set  is  a  collection  of  points  such  that,  if  x  and  y 
are  any  two  points  in  the  collection,  every  point  on  the  line  segment 


15 


joining  x  and  y  is  olso  in  the  collection.  Any  point  in  a  convex  set 
which  does  not  lie  on  the  line  segment  joining  some  two  other  points 
in  the  set  is  known  as  an  extreme  point.   (See  Figure  4.)  .  A  rectangle, 
for  instance,  is  a  convex  set  and  the  four  corners  are  extreme  points. 
A  convex  set  which  has  a  finite  number  of  extreme  points  is  given  the 
special  name  of  convex  polyhedron. 


Extreme 
Point 


Convex  Set  Non-Convex  Set 


Fig.  U     Convex  and  Non-Convex  Sets 

In  Figure  5  the  convex  polyhedron  P.  P  p  p  is  said  to  gen- 
erate the  cone  0-P,  P?  p..  P..  This  means  that  those  rays  emanating 
from  the  origin  0  which  pass  through  the  perimeter  of  the  convex 
polyhedron  P,  P^  P_  P.  form  the  sides  of  a  cone.  Such  a  cone  is  re- 
ferred to  as  a  convex  polyhedral  cone . 
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x2 


Fig.  5  Convex  Polyhedral  Cone 

Clearly  the  unit  vectors  in  an  m-dimensional  space  generate  a 
convex  polyhedral  cone  vhich  will  encompass  all  positive  vectors  in 
the  m-d linens ional  space. 

The  linear  transformation  of  points  has  been  discussed.  Just 
as  points  in  an  n-dimensional  space  can  be  transformed  into  image 
points  in  an  m-dimensional  space,  so  can  a  convex  polyhedron.  Indeed, 
the  image  under  (resulting  from)  a  linear  transformation  of  a  convex 
polyhedron  is  a  convex  polyhedron.  An  illustration  of  this  is  shown 
in  Figure  3 .  The  convex  polyhedron  in  space  N  formed  by  the  points 

e  "  t  e^'),   and  e^'  (a  triangle)  has  been  transformed  into  a  triangle 
formed  by  the  points  y(  ',  y(?)y  and  y(3)  in  space  M 


Simplex  Method  of  Solution 
The  preceding  mathematical  discussion  has  presented  some  basic 
information  necessary  to  the  understanding  of  what  a  linear  programming 
problem  is  and  how  it  can  be  solved.  A  simple  numerical  example  will 
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now  be  presented  and  its  solution  will  be  obtained  by  the  simplex 
method. 

The  example  consists  of  two  inequations  in  three  unknowns  and 
a  linear  function  (sometimes  called  a  linear  functional)  of  the  three 
unknowns.  The  linear  functional  is  to  be  maximized.  A  restriction  is 
imposed  that  the  unknowns  when  found  must  be  non-negative,  that  is, 
must  equal  or  be  greater  than  zero  (20). 

2x4-    y  *  2  z  &  JO 

*  ♦  2  y  ♦  3  z  £.J4 
f  (*>  y»  z)   =  x  <*•  y  +  z   (maximum) 
with  x^O;  y>0j  z^O 
To  convert  the  inequalities  to  equations  a  non-negative  variable 
is  added  to  each  ineouation  and  the  model  appears  as  follows: 

u         4.2x       +       y4.2zal0 
1. 

v     4-  x       +2y       t       3z     =     U 

2.  x,   y,   z,   u,   v,    >    0 

3.  f  (x,  y,  z)  s  x  +  y  4-  z  4-   (0)  u  *  (0)  v  to  be  a  maximum. 
Stated  in  words  the  problem  is  to  find  the  point  or  points  in  a 

5-dimensional  space  U  which  have  non-negative  coordinates  (u,  v,  x,  y,z) 
satisfying  equations  (1)  above  and  such  that  the  linear  functional 
(x  +  y  +  z)  is  a  maximum.  It  can  be  shown  that  the  point  or  set  of 
points  so  described  forms  a  convex  set.  The  coefficients  associated 
with  the  variables  in  equations  (1)  plus  the  right  side  of  the  equa- 
tions yield  6  column  vectors  P.  in  the  2-dimensional  space  W.  Thus, 
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pi-'o)>  >,.#),  P3=f12)s  ^-(i), »,.(!),  andP0.;a 

Equations  (1)  can  be  restated: 

Pl        P2       F3     F4       p5     p0 

•(c)  ♦  *'?;  t  *[!)  r{\)  ♦ .(«)  .  g] 

or  uPx  +  vP2  |  xP3  4-  yP^  +  ZP5  -  po 

Recall  the  discussion  of  linear  transformations.  The  coefficients 
in  the  above  equations  form  a  matrix  associated  with  the  linear  trans- 
formation which  takes  the  point  Q  ■  (u,  v,  x,  y,  z)  in  5-dimensional 
space,  U,  into  the  point  Fn  =  (10,  H)  in  ^-dimensional  space,  V. 

In  the  equations  above,  u  is  the  coordinate  along  axis  1,  v 
along  axis  2,  in  the  space  U. 

T  .<«  =  Pj  .  (1)  g<«  ♦  (0)  g(2> 

T  e(:>)  5  P2  .  (0)  gU)  ♦  (1)  g(2) 

T  e<3>  „  P3  .  (2)  g««  *  (1)  g(« 

T  e(i>  =  P4  .  (1)  g(l)  +  (2)  g(2) 

T  e(5)  =  P.  =  (2)  gd)  +  (3)  g(2)  vhere  the 

coefficients  of  the  g's  are  the  elements  in  the  re- 
spective vectors  P.. 

There  is  more  than  one  point  Q  in  U  having  non-negative  co- 
ordinates u,  v,  x,  y,  z,  which  will  satisfy  the  equations,  i.e., 
satisfy  T  (Q)  ■  Pq,  and  it  can  be  shown  that  these  values  of  Q   form 
a  convex  set  in  U.  For  example:  Q  =  (u,  v,  x,  y,  z)  = 

15.  Ibid.  p.  53. 
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(2/3,  0,  0,  0,  H/3);  Q  ■  (10,  U,  0,  0,  0);  and  Q  =  (0,  C,  1/2,  0,  9/2), 
are  each  answers  which  will  satisfy  the  two  equations,  as  substitution 
back  into  the  equations  will  show.  All  possible  Q'a  which  have  non- 
negative  coordinates  u,  v,  x,  y,  and  z  form  a  convex  set  which,  because 
it  has  a  finite  number  of  extreme  points  is  a  convex  polyhedron. 

An  important  theorem  states  that  a  linear  functional  f  (u,  v,  x, 
y,  z)  defined  on  a  convex  polyhedron  takes  on  its  maximum  (or  minimum) 
at  an  extreme  point  of  the  convex  set.    If  the  extreme  points  of  the 
convex  polyhedron  generated  by  all  values  of  Q  which  have  non-negative 
coordinates  could  be  found,  the  maximum  value  of  the  functional  could 
also  be  found.  That  is  to  say,  having  the  coordinates  of  an  extreme 
point,  and  substituting  them  into  the  functional,  a  value  of  the  func- 
tional would  be  obtained.  Comparing  this  value  of  the  functional  with 
the  value  obtained  from  every  other  extreme  point  of  the  convex  set, 
one  could  pick  out  the  extreme  point  whose  coordinates  yielded  the  max- 
imum value  of  the  functional.  These  coordinates,  or  values  for  x,  y, 
and  z  (u  and  v  are  of  no  interest,  having  been  added  merely  to  chanre 
the  inecuations  to  equations)  are  then  the  required  values  for  the 
optimum  solution  of  the  problem. 

The  procedure,  then,  is  to  locate  an  extreme  point,  calculate 
the  value  of  its  functional,  proceed  to  another  extreme  point,  and  so 
on,  until  the  reouired  maximum  value  of  the  functional  is  obtained. 
The  extreme  point  which  gives  the  maximum  value  of  the  functional  is 
the  solution. 

16.  Ibid.  P.  52. 
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Recalling  the  discussion  of  linear  independence,  it  was  stated 
that  in  2-dimensional  space  the  maximum  number  of  points  or  vectors 
in  a  linearly  independent  set,  or  basis,  is  two.  The  above  points  P 
are  all  in  2-cimensional  space  and  subjecting  any  pair  of  points,  P,  to 
the  criterion  for  linearly  independent  points  as  described  on  page  &. 
it  can  be  seen  that  any  combination  of  two  such  points  P  is  a  linearly 
independent  set.  Any  point,  therefore,  can  be  expressed  as  a  linear 
combination  of  any  two  other  such  points  P.  For  example,  if  x,  y,  and 
z  were  each  eoual  to  zero,  the  ecuction  could  be  stated  as 
u  ?1     +  v  ?2     +  OP.  +  OP   4-  0  P5  =  PQ 

or  u  P1  +  v  ?2     s  P0,   Inserting  the  coordinates  of 
?lt     P2,  and  PQ 

u  (1)   ♦  v  (0)   «  10  or  u  .  10 
u  (C)  4.  v  (1)  -  Uor  v  ■  H 
The  point  Q  ■  qq  -j,'  ?'  ?'   ^\  is  an  extreme  point.  This 

point  Q  fulfills  the  requirements  of  the  equations  of  restrictions,  but 
the  value  of  the  functional  resulting  from  this  point  C  is  not  a  max- 
imum, but  is,  instead,  zero.  The  simplex  method  is  designed  to  permit 
movement  from  one  extreme  point  to  another  and  in  such  a  direction  that 
an  increase  in  the  functional  will  result. 

To  demonstrate  the  simplex  method  for  the  solution  of  this 
problem  the  P  's  are  placed  in  a  box  arrangement.  The  P~  has  been 
moved  into  the  first  position. 
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Table  1 


po 

pl 

P2 

P3 

h 

P5 

10 

1 

0 

2 

1 

2 

H 

0 

1 

1 

2 

3 

The  original  equations  of  the  problem  can  be  derived  from  the 
box  arrangement  as  follows:  The  single  vertical  lines  represent  plus 
signs  and  the  double  lines  separating  PQ  from  the  rest  of  the  arrange- 
ment represent  equal  signs.  Visualizing  the  elements  of  eech  horizon- 
tal row  as  being  multiplied  by  their  respective  unknowns  and  including 
the  equal  and  plus  signs  one  can  recreate  the  two  equations.  Thus  from 
row  1: 

10  =  u(l)   f  v(C)   ♦  x(2)   *  y(l)   *  z(2) 
and  from  row  2: 

H  =  u(0)   4.  v(l)   4.  x(l)   4.  y(2)   *  Z(3) 

An  initial  solution  was  obtained  utilizing  the  basis  set  P-^,  P2 
To  show  this,  Table  1  is  expanded  into  Table  2,  where  P^  and  P2  have 
been  placed  to  the  left  of  Pq  in  the  stub  column. 

Table  2 


c  -*- 

1 

0 

0 

1 

1 

1 

p5 

\ 

Vector 

P0 

pl 

p2 

p3 

h 

0 

Pl 

10 

1 

0 

2 

1 

2 

0 

P2 

u 

0 

1 

1 

2 

3 

Table  2  has  the  property  of  expressing  any  of  the  vectors  in 
terms  of  the  set  P-^  and  P2«  Any  vector  as  designated  by  the  column 
headings  can  be  stated  by  multiplying  each  element  in  the  column  by 
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the  vector  on  its  row  in  the  stub  find  adding.  Thus  for  Pq  the  element 
on  the  first  line  under  PQ,  10,  is  multiplied  by  P1  and  added  to  the 
product  of  the  second  element,  14,  with  Pp. 
Thus  P0  =  (10)  ?1     +  (U)  P2 

Similarly   Px  .  (1)  ?1     *  (0)  P2 

V2     «  (0)  ?1  +  (1)  P2 
P3  «  (2)  P-l  *  (1)  P2 

P^  .  (1)  P1  *  (2)  P2 

P5  «  (2)  PX  ♦  (3)  P2 

This  follows  from  the  fact  that  any  point  in  a  2 -dimensional 
space  can  be  expressed  as  a  linear  combination  of  the  points  in  a 
linearly  independent  set.  Any  of  the  above  can  be  verified  by  inserting 
the  numerical  values  of  the  P^'s.  Thus  for  P, 

[l]  -   m  (5)    *<2>(?) 

The  numbers  in  the  top  row  and  far  left  column  designated  by  c 
are  the  coefficients  of  the  corresponding  variables  in  the  functional 
whose  maximum  is  sought.  Thus  the  number,  1,  above  P-,  P  ,  and  P^.  is 
the  coefficient  in  the  functional  of  x,  y,  and  z,  respectively.  The 
coefficient  of  u  and  v  is  0  as  shown  by  the  0  above  and  to  the  left  of 
Px  and  P2. 

Table  2  completes  the  first  set  of  calculations.  All  of  the  six 
vectors  (P^,  P2,  P^>  P,»  P*,  and  Pq)  appearing  at  the  top  of  the  table 
have  been  solved  or  stated  in  terras  of  the  basis  P^,  P2  and  a  feasible 
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solution  is  indicated  by  the  elements  under  Pnj  ie.  u  =  10,  v  ■  H. 
Because  P-,  P.,  and  P<-  do  not  appear  in  the  stub  their  corresponding 
variables  x,  y,  and  z  are  each  equal  to  zero.  The  value  of  the  func- 
tional at  this  stage  is  obtained  by  multiplying  the  feasible  solution 
listed  under  P0  by  the  functional  coefficients  listed  to  the  left  under 
c.  The  value  of  the  functional  is  thus  (10)  (0)  4-  (H)(0)  -  0. 

As  was  stated,  it  is  desired  to  move  to  another  extreme  point  to 
investigate  the  value  of  the  functional.  The  simplex  method  provides 
an  algorism  which  yields  a  new  value  of  the  functional  that  is  no  smal- 
ler  than  the  previous  value.  That  is  to  say,  the  new  value  of  the  ^ 
functional  will  either  be  larger  than  the  previous  value  or  possibly 
the  same  value  as  the  previous  value  of  the  functional  but  will  never 
be  smaller.  The  new  extreme  point  found  will  increase  the  value  of  one 
of  the  variables  x,  y,  or  z  at  the  expense  of  decreasing  the  present 
values  of  u  and  v  namely  10  and  14-,  respectively.  To  find  which  vari- 
able to  increase  (what  vector  Po,  P/,  or  ?c   to  bring  into  the  basis), 
the  unit  functional  weights  of  x,  y,  and  z  are  computed  and  compared. 
The  largest  value  will  designate  the  vector  to  bring  into  the  basis. 

For  example,  an  increase  in  the  size  of  x  by  one  unit  (from  0  to 
1)  will  increase  the  size  of  the  functional  by  1,  the  functional  coef- 
ficient of  x.  But  from  Table  2,  1  P3  3  2  P-^  +  1  P2  or  an  increase 
in  x  by  1  will  decrease  u  by  2  and  v  by  1.  The  decrease  in  the  func- 
tional caused  by  reduction  of  u  and  v  is  equal  to  (0)(2)  *  (0) (1)  =  0 
because  the  functional  coefficients  of  u  and  v  are  each  zero.  A  unit 
increase  of  x,  therefore,  will  cause  the  functional  to  increase  in 
value  by  1. 
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Similarly,  a  unit  increase  of  y  will  cause  the  functional  value 
to  increase  by 

cy   y   Cu   u   c^   v 

(1)  (1)  -  (0)  (1)  -  (0)  (2)  »  1 

The  remaining  variable  z,  if  increased  by  1,  will  likewise 
increase  the  functional  by  1. 


cz   z 


cv   V 


Cu   u 
(1)  (1)  -  (0)  (2)  -  (C)  (3)  »  1 

Since  the  unit  functional  weights  of  x,  y,  and  z  are  equal,  the 
value  of  the  functional  will  be  increased  the  same  amount  no  matter 
which  variable  is  increased  (by  inclusion  of  its  F  in  the  stub  basis) , 
hence  it  is  immaterial  which  one  is  chosen.  If  the  unit  functional 
weights  are  not  eoual,  the  one  with  the  largest  positive  value  desig- 
nate s  the  F  to  include  in  the  basis. 

Table  2  provides  a  means  for  a  systematic  computation  of  the 
unit  functional  weights  as  will  now  be  shown.  The  table  is  reproduced 
here  for  convenience  as  Table  3. 

Table  3 


Cnv 

1 

0 

0 

1 

1 

1 

1 

Vec- 
tor 

P0 

pl 

p2 

F3 

h 

P5 

0 

Pl 

10 

1 

0 

2 

1 

2 

0 

p? 

u 

0 

1 

1 

2 

3 

F 

0 

0 

0 

0 

0 

0 

C-F 

0 

0 

1 

1 

1 

The  values  in  the  row  designated  by  F  are  obtained  as  follows: 
each  element  in  a  column  is  multiplied  by  the  c  which  appears  to  the 
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far  left  in  the  element's  row  and  the  products  are  added  to  give  the 
value  of  F  for  that  column.  Thus  the  value  of  F  under  Fn  Is  equal  to 
(10)  (0)   4-  (H)(0)   =  Oj  F  under  P3  is  equal  to  (2)(0)  ±  (1)  (0)  -  0, 
etc.  The  F  value  under  PQ  is  the  value  of  the  functional  obtained  from 
the  solution  Q  -  (10,  H,  0,  ?,  0)  .  The  F  value  in  each  of  the  other 
columns  is  the  value  of  the  decrease  of  the  functional  which  would  re- 
sult from  the  decrease  in  value  of  the  two  variables  u  and  v  if  a  unit 
increase  in  the  variable  corresponding  to  that  column  were  made.  Thus, 
if  x  were  increased  from  0  to  1,  the  value  of  F_  =  0  indicates  that  the 
reductions  of  u  and  v  would  have  no  effect  upon  the  value  of  the  func- 
tional. 

The  elements  in  the  row  C-F  are  the  unit  functional  weights,  and 
are  obtained  by  subtracting  the  F  value  in  each  column  from  the  corres- 
ponding C  value  at  the  top  of  the  column. 

The  unit  functional  weights  (C-F)  are  examined.  The  largest 
positive  value  determines  which  vector  F  to  introduce  into  the  basis  in 
the  stub  to  replace  one  P  which  is  now  in  the  stub.  This  replacement 
of  an  existing  P  in  the  basis  by  a  new  one  corresponds  to  moving  from 
one  extreme  point  to  another.  By  introducing  the  P  which  has  the 
largest  unit  functional  weight  (C-F)  the  largest  increase  in  the  value 
of  the  functional  is  insured.  As  long  as  there  are  positive  values  of 
(C-F)  an  increase  in  the  functional  is  possible. 

In  the  present  problem  there  are  three  equal  positive  values  of 
(C-F),  each  equal  to  1.  It  is  immaterial  v/hich  one  is  chosen  though 
the  length  of  the  computations  might  be  shortened  somewhat  by  a  judi- 
cious choice.   For  illustration,  P^  is  chosen  to  replace  one  of  the 
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vectors,  P^  or  P^,  in  the  stub. 

The  choice  of  which  vector  P1  or  P^  in  the  stub  to  be  replaced 

by  the  incoming  vector  P^  is  made  as  follows:  Divide  each  positive 

element  in  the  Tr   column  into  the  corresponding  element  in  the  Fq  column 

and  compare  the  results .  Thus 

10     , 
~2    =  5 

lj    s  U  1/3 

The  smallest  result,  here  U   1/3,  designates  the  vector  P2  of  the 
stub  as  the  one  to  be  replaced  by  Pr.  U   1/3  or  (14/3)  is  the  new  value 
of  z.  Since  the  vector  P0  is  being  replaced,  its  variable  v  decreases 
in  value  from  14  to  C.  A  blank  table  similar  to  Table  3  is  now  lai^ 
out  but  with  ?2   *n  the  stub  replaced  by  Pc.  Note  that  the  c  •  1  0f  Pe- 
ls placed  to  the  left  of  F5  in  the  stub  of  Table  A. 


Table  U 


c> 

0 

0 

1 

1 

1 

1 

Vec- 
tor 

P0 

Pl 

P2 

P3 

?A 

0 

pi 

1 

p5 

To  complete  the  table,  one  recalls  that  a  property  of  this  ar- 
rangement is  that  any  vector  as  listed  at  the  head  of  a  column  can  be 
expressed  as  a  linear  combination  of  the  two  vectors  in  the  stub.  If 
the  vectors  P-^,  P2,  P3,  P^,  P*,  can  be  expressed  as  linear  combinations 
of  P]_  and  ?ct   the  coefficients  will  be  the  elements  of  the  table.  From 
Table  3  the  following  equations  hold: 
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(a) 

po 

r 

10  ?1 

X 

(b) 

pl 

= 

1?1 

4> 

(c) 

p2 

r 

0P1 

4 

(d) 

P3 

a 

2P1 

t 

(e) 

P4 

5 

1P1 

4- 

(f) 

P5 

: 

2P1 

+ 

0  P2 


1  F« 


1  P- 


2  P. 


3  Pi 


From  equation  (f)  P2  =  -2/3  Px  4-  1/3  ?5   thus  -2/3  and  1/3  are 
the  elements  in  the  new  table  under  P„.  Substituting  -2/3  ?i   *  1/3  p^ 

into  equation  (a) 

P0  =  10  ?1    f  H  (-2/3  Px  4-  1/3  P5) 

or       P0     =     2/3  ?i  4-     M  p..     Thus   2/3  and  U/3  are  the  elements 
to  be  filled  in  under  P   . 


Similarly, 
P 


lpi  + 


0  F, 


1     -  *  *1     x         w  *5 

P3    =      4/3  Px    ♦    1/3  P5 

P^     =    -1/3  P1     4-     2/3  P5 


0P1     + 


1  Pi 


determine  the  elements  for 


the  other  columns.     The  values  of  F  and  C-F  are  determined  in  the  same 
manner  as  employed   in  Table  3.     The  completed  table,   Table   5  then 
appears  as   follows : 
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Table  5 


0^- 

1 

0 

0 

1 

1 

1 

♦ 

Vec- 
tor 

P0 

pl 

?2 

P3 

P4 

P5 

0 

Pl 

2/3 

1 

-2/3 

V3 

-1/3 

0 

1 

P5 

U/3 

0 

1/3 

1/3 

2/3 

1 

F 

U/3 

0 

1/3 

1/3 

2/3 

1 

C  -  F 

0 

-1/3 

2/3 

1/3 

0 

Table  5  shows  (in  the  PQ  column)  that  a  new  solution  u  =  2/3, 
z  =  L+/3,  with  x=0,  y"0,  v"0,  has  been  obtained.  The  F  value 
below  P  shows  that  the  value  of  the  functional  has  been  increased 
from  0  to  m/3,  and  the  two  positive  values  in  the  C-F  row  show  that 
further  improvement  is  possible.  The  2/3  value  of  C-F  in  column  P3 
shows  that  P^  is  the  most  profitable  vector  to  introduce  into  the 
basis.  As  before  the  division  of  the  elements  under  P.  by  those 
under  P~  shows  that  Pj  is  the  vector  in  the  besis  to  be  replaced  by 
P3.  Carrying  out  the  techniques  employed  previously  Table  6  is 
constructed  using  the  equations  obtained  from  Table  5  to  compute  the 
elements  in  Table  6. 
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Table  6 


e-*- 

0 

0 

1 

1 

1 

T 

Vec- 
tor 

P0 

Pl 

P2 

P3 

*u 

p5 

1 

P3 

1/2 

3/4 

-1/2 

1 

-iA 

0 

1 

P5 

9/2 

-iA 

1/2 

0 

3A 

1 

F 

5 

1/2 

0 

1 

1/2 

1 

C  -  F 

-1/2 

0 

0 

1/2 

0 

Examination  of  Table  6  reveals  a  new  feasible  solution  for  Q, 
namely,  Q  »  (0,  0,  l/2,  0,  9/2),  and  an  increased  value  of  the  function- 
al, namely  5.  The  positive  value  of  C-F,  the  l/2  in  column  P.,  indi- 
cates that  a  larger  value  of  the  functional  can  be  obtained  by  introdu- 
cing P^  into  the  basis.  Because  3A  in  row  P^  is  the  only  positive 
value  in  the  P^  column,  P^  is  the  vector  in  the  basis  to  be  replaced. 
Carrying  through  the  same  procedure  as  before  by  computing  the  elements 
for  the  next  table  from  the  ecuations  indicated  in  Table  6,  the  final 
table,  Table  7  is  completed 

Table  7 


( 

*  > 

0 

0 

1 

1 

1 

* 

' 

Vec- 
tor 

p0 

pl 

P2 

p3 

h 

P5 

1 

P3 

2 

2/3 

-1/3 

1 

0 

1/3 

1 

h 

6 

-1/3 

2/3 

0 

1 

V3 

F 

8 

1/3 

1/3 

1 

1 

5/3 

C-F 

-1/3 

-1/3 

0 

0 

-2/3 
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The  non-appearance  of  any  positive  values  of  C-F,  the  unit 
functional  weights,  indicates  that  no  further  increase  in  the  value 
of  the  functional  is  possible.  Reading  in  the  Pn  column  it  is  seen 
that  the  optimum  solution  is  x  =  2,  y  ■  6,  u  «  0,  v  ■  0,  z  B  0  and  the 
maximum  value  possible  of  the  functional  is  8. 

The  procedure  for  arriving  at  the  table  elements,  namely  by 
solution  of  equations  from  the  preceding  table,  becomes  extremely  labor- 
ious when  the  number  of  variables  and  the  number  of  equations  is  in- 
creased. A  method  exists  that  permits  easy  calculation  of  the  elements 
by  a  systematic  mechanical  method.  This  method  permits  the  solution  of 
extremely  lerge  problems  involving  many  variables  and  many  equations  in 
the  original  problem. 

The  sample  problem  which  has  been  solved  originated  as  a  set  of 
two  inequations.  3y  introducing  the  additional  unknowns  to  convert  the 
inequations  to  equations  the  vectors  P-^  c(n)  *nd  Pp  z[i)   were  intro- 
duced. These  two  slack  vectors,  as  they  are  called,  provided  a  known 
basis  from  which  the  first  feasible  solution  (u  =  10,  v  =  14)  was  ob- 
tained. When  the  si2e  of  the  problem  is  increased  the  introduction  of 
such  a  simple  basis  becomes  of  more  importance.   If  the  original  set  of 
restrictive  conditions  are  equations  instead  of  inequations,  slack 
vectors,  now  known  as  artificial  vectors,  may  still  be  added  to  provide 
a  basis.  The  associated  variables,  however,  must  be  included  in  the 
functional  to  be  maximized  and  each  such  variable  must  have  a  negative 
coefficient  of  sufficient  size  so  as  to  insure  that  the  variable  will 
not  appear  in  the  final  answer. 

The  above  sample  problem  was  concerned  with  the  maximization  of 
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a  linear  functional.  As  was  stated  previously  the  optimization  in  a 
linear  programming  problem  may  be  the  minimization  of  some  function 
such  as  cost.  A  minimization  problem  may  be  solved  in  a  similar  fashion 
except  that  successive  extreme  points  are  found  by  an  algorism  vhich 
insures  that  each  new  value  of  the  functional  thst  is  obtained  will  be 
no  larger  thon  any  previous  one. 
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ASSUMPTIONS,  POSSIBILITIES  AND  LIMITATIONS 
The  bibliography  contains  accounts  of  the  application  of  linear 
programming  to  industrial  problems.  In  the  main,  however,  the  practi- 
cality of  the  technique  suffers  from  certain  restrictive  assumptions, 
difficulties  of  model  building,  lack  of  proper  data,  necessity  for 
complex  computing  machines  and  the  necessity  that  a  linear  programming 
program  be  set  up,  if  not  operated,  by  highly  trained  personnel.  The 
last  two  may  impose  severe  financial  restrictions. 

The  assumptions  that  must  be  met  in  order  to  apply  linear  program- 
ming are  as  follows:  ' 


1.  The  assumption  of  linearity.  A  linear  process  is  one  whose  s 

18 

ratio  of  output  to  input  is  constant.    Thus  for  an  operating  machine 

tool,  if  an  input  of  one  unit  of  raw  material  produces  one  unit  of 
output  in  a  unit  of  time,  an  input  of  100  units  must  produce  100  units 
of  output  in  100  units  of  time .  Such  en  assumption  neglects  such 
aspects  as  rejects,  tool  wear,  delay  in  delivery  of  materia],  set-up 
time,  industrial  fatigue,  or  unauthorized  absence  of  the  operator,  etc. 
Such  non-linearity  of  a  production  process  is,  of  course,  nothing  new 
and  confronts  production  control  personnel  no  matter  what  method  they 
employ.  Because  of  the  large  number  of  computations  made  in  linear 
programming,  however,  and  the  fact  that  the  results  or  decisions  may 


17.  Dorfman,  R.,  Application  of  Linear  Programming  to  the  The ory  of  the 
Firm;  Berkeley,  University  of  California  Press,  1951,  p.  80. 

18.  This  ratio  is  a   result  of  the  teclinical  aspects  of  the  specific 
process  or  activity.  In  the  machining  of  metal,  for  instance, 
specified  speeds  and  feeds  in  addition  to  other  specifications 
determine  the  process.  If  the  speed  or  feed  were  changed  a 
technically  feasible  process  might  result.  Such  a  modification  of 
a  process, however,  is,  in  linear  programming  terms  ,  a  new  and 
distinct  process. 


? 
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affect  a  large  segment  of  a  firm's  production  program  it  becomes  neces- 
sary to  insure  a  high  degree  of  accuracy.  Morgenstern  states: 

Linear  programming,  or  any  other  similar  utilization  of  great  masses 
of  econom.ic  data,  cannot  be  expected  to  make  decisive  practical 
progress  until  there  is  satisfaction  that  the  data  warrant  the  im- 
plied extensive  and  costly  numerical  operations .19 

Other  examples  of  non-linearity  are  the  reduction  of  sales  price 

per  unit  for  volume  sales  and  the  rise  in  productivity  resulting  from 

minor  technological  improvements. 

2.  The  assumption  of  divisibility.  Any  process  may  be  operated 
at  any  positive  level.  This  assumes,  for  instance,  that  if  one  large 
press  can  stamp  out  100  fenders  in  an  hour  then  one  half  of  a  press  can 
stamp  out  50  fenders  in  an  hour. 

3.  The  assumption  of  additivity: 

It  is  assumed  that  two  or  more  processes  can  be  used  simultaneously, 
within  the  limitations  of  available  resources,  and  that  if  this  is 
done  the  quantities  of  the  outputs  and  inputs  will  be  the  sums  of  the 
quantities  which  would  result  if  the  several  processes  were  used 
individually.20 

L>.   The  assumption  that  the  number  of  processes  is  finite.   Be- 
cause of  the  restricted  field  of  interest  of  many  firms  the  processes 
or  activities  are  indeed  finite.  Some  industries,  however,  such  as 
agriculture,  oil  refining,  or  the  chemical  industries  have  an  infinite 
range  of  choices . 

Melvin  E.  Salveson  has  listed  some  of  the  problems  which  can  and 


19.  Morgenstern,  0.,  "The  Accuracy  of  Economic  Observations,"  Activity 
Analysis  of  i'rodurtlon  and  Allocation,  Cowles  Comm.  for  Research  in 
Economics,  Monograph  13;  New  York,  John  Wiley  and  Sons,  1951,  p.  283. 

20.  Dorfman,  R.,  Application  of  Linear  Programming  to  the  Theory  of  the 
Firm;  Berkeley,  University  of  California  Press,  1951,  p.  81. 
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21 
cannot  be  solved  by  linear  programming.    He  states  that  this  tech- 
nique can  be  used: 

1.  to  determine  the  optimum  shop  load  during  any  one  time  period  or 
over  several  time  periods. 

2.  to  determine  the  optimum  amount  of  overtime  to  use  on  any  item  or 
machine  tool. 

3.  to  determine  the  optimum  mix  of  commodities  to  make  in  the  shop  or 
factory. 

4.  to  determine  the  amount  of  rerouting  of  work  in  the  shop. 

5.  to  determine  the  optimum  level  of  inventory  (raw,  in-process,  and 
finished) . 

6.  to  determine  the  optimum  distribution  of  production  of  commodities 
between  time  periods. 

He  states  that  linear  programming  cannot  be  used: 

1.  to  determine  a  schedule  for  a  shop. 

2.  to  select  the  best  sequence  of  production  in  order  to  minimize, 
say,  set  up  time. 

3.  to  eliminate  need  for  expediters. 

U.   to  avoid  all  conflicts  in  production,  such  as  temporary  bottle- 
necks and  hence,  its  program  is  not  necessarily  always  achievable 
(unless  it  is  constructed  with  overly  liberal  delay  allowances) . 

5.  to  give  delivery  dates  more  precisely  than  the  length  of  its  time 
periods. 

The  reasons  for  the  non-applicability  of  linear  programming  to 
the  five  items  above  are  as  follows: 

(a)  Item  1  and  item  5.  These  refer  to  the  scheduling  of  activi- 
ties. Production  or  any  other  activity  of  a  linear  programming  problem 
when  considered  in  relation  to  time  is  considered  to  begin  at  the  start 
of  a  specified  time  period,  continue  during  the  time  period,  at  a  con- 
stant rate  and  to  reach  completion  at  the  end  of  the  time  period.  This 
activity  time  requirement  is  faced  in  machine  or  shop  loading  problems 
also  but  suitable  assumptions  or  restrictions  in  the  formulation  of  the 


21.  Salveson,  M.E.,  "Mathematical  Models  in  Management  Programming," 
Journal  of  Industrial  Engineering.  March  1954,  vol.  5,  No.  2, 
p.  10. 
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mathematical  model  may  allow  an  operationally  meaningful  solution. 
An  example  of  such  a  restrictive  assumption  will  be  given  later. 

(b)  Item  2.  The  selection  of  the  sequence  of  production  is 
determined  by  the  technology  of  the  production  processes.  Thus,  for 
example,  the  manufacture  of  a  certain  item  might  require  that  a  milling 
operation  precede  a  drilling  operation.  This  seauence  is  part  of  the 
technological  requirements  for  the  manufacture  of  the  item.  In  addi- 
tion ,  set  up  time  is  a  non-linear  activity  or  function  of  a  process. 
Thus  one  hour  of  set  up  time  may  be  required  for  an  operation  regard- 
less of  the  level  of  activity  of  the  process. 

(c)  Item  3  and  item  4.  These  restrictions  on  the  applicability 
of  linear  programming  to  production  result  from  the  inability  of  any 
mathematical  model  to  account  fully  for  the  emergencies  of  the  every 
day  world.  If  every  exigency  were  known  ahead  of  time  it  might  pos- 
sibly be  incorporated  in  the  model.  The  possibility  of  forecasting 
the  occurence  of  every  conceivable  perturbation  or  disturbance  in  the 
orderly  activities  of  nature  is,  of  course,  extremely  slim.  Many 
conflicts  in  production,  such  as  machine  breakdown  or  non-delivery  of 
raw  material  due  to  personnel  error,  may  be  compensated  for  by  con- 
structing the  model  with  very  liberal  delay  allowances.  Expediters 
are  necessary  to  keep  the  actual  production  as  close  to  the  program  as 
possible. 
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MODELS 

Many  mathematical  models  heve  been  presenter?  in  linear  program- 
ming literature  and  work  in  the  area  of  model  building  is  proceeding. 
To  show  some  of  the  range  of  possibilities  a  few  models  are  presented 
here. 

Machine  Assignment^ 

A  number  oi  operations  are  to  be  performed  in  each  of  a  number 
of  mutually  exclusive  operation  types.  A  number  of  machines  of  various 
types  are  available.  The  productivity  of  each  machine  type  in  each  of 
the  various  types  of  operations  is  known . 

Required:  to  assign  the  machines  to  operations  such  that  the 

total  productivity  of  all  machines  on  all  operations  is  a  maximum. 

Let: 

m  =  the  number  of  types  of  machines. 

a^  =  the  number  of  machines  of  type  (i)  where 
i  =  1,  2,  3}  •  •  i  •)  m. 
n  ■»  the  number  of  types  of  operations, 
b.  =  the  number  of  operations  of  type  (j)  where  j  = 

1>  2,  3, ,  n. 

Cji   =  the  productivity  of  a  machine  of  type    (i)    in 

performing  operation   (j). 
xi1  =  ^e  nUFj^er  °^  type   (i)   machines  to  be  assigned  to 
a  type   (j)    operation. 


22.  Harrison,  J.O.,  Jr.,  Linear  Programming  and  Operations  Research, 
Informal  Seminar  in  Operations  Research,  Paper  No.  2;  Baltimore, 
Johns  Hopkins  Univ.,  October  13,   1953,   pp.   12-13 • 
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The  x.j  are  to  be  selected  such  that: 


£    0 

n 

3*/ 

m 

< 

ai 

>     *U 

- 

b0 

1='  v*^ 

Production  Line  Loading  to  Meet  Sales  Requirements  ^ 
Given:  one  production  line  producing  one  type  of  product  which 
sells  for  a  fixed  unit  price.  The  unit  costs  of  regular  time  produc- 
tion, overtime  production,  and  storage  are  known.  The  rates  of  produc- 
tion per  unit  time  are  known.  A  sales  demand  during  each  of  a  number 
of  successive  time  periods  is  known.  The  time  period  is  one  month. 

Desired:  a  production  loading  plan  which  will  meet  the  sales 
demand  and  minimize  the  combined  costs  of  production  and  storage. 
Inventory  is  taken  at  the  end  of  each  month. 

Let: 

k   s  number  of  time  periods  (months)  to  load. 

SJ  r  number  of  units  of  finished  product  to  be  sold  during  the 

ith  period. 
Iq  -  initial  inventory. 
R  z  maximum  number  of  units  which  can  be  produced  during  a 

month  on  regular  time. 
U   =  maximum  number  of  units  which  can  be  produced  during  a 

month  on  overtime. 
CT  -  cost  of  storage  of  one  unit  of  product  for  one  month. 


23.  Ibid.  pp.  M-15. 
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Cx  =  unit  production  cost  utilizing  regular  time. 
Cy  -  unit  production  cost  utilizing  overtime. 
Find:  . 

x^   s  number  of  units  to  be  produced  on  regular  time  during 

the  ith  month. 
y^  r  number  of  units  to  be  produced  on  overtime  during  the 

ith  month. 

(a)  x±  ,  y±  Z    0         (i  .  1,  2, ,k) 

(b)  xi  <   R 

(c)  y-L  <  U 

i  i 

(d)  ^(xa  +  ya)   >   V^sa   -I0 

k  ' 

(e)  \   (Cj  I±     4.  Cx  x^  f  C_  y^)  to  be  minimized,  where 

and  (a  =  1,  2, ,  i) 

Equation  (d)  is  obtained  as  follows: 

The  inventory  during  any  month,  a,  is  denoted  by  I&.  Since 
the  inventory  at  the  end  of  any  month  is  equal  to  the  inventory  at  the 
end  of  the  previous  month  plus  the  monthly  production  minus  the  monthly 

sales, 

xa  ■  xa-l  +  xa  *  ya  ~  sa 

Summing  this  equation  from  a  *  1  to  a  =  i  yields  equation   (f). 

Since  the  right  side  of  equation   (f)   must  be  equal  to  or  greater  than 

zero  (  >  0) ,  or  j.  i 

I0     +       V(xa     +     ya)      -    Y%a     >       0 

A  slight  rearrangement  yields  equation  (d) . 
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A  Numerical  Example  -  Product 
Assembly  Over  Consecutive  Time  Periods'"* 

To  illustrate  the  complexity  to  which  even  problems  of  small 
size  are  subject  a  numerical  example  is  presented. 

Assumptions:  a  factory  making  three  products,  each  of  which 
can  be  final  products  and  two  of  which  are  in  addition  intermediate 
products . 

Product  1  (F^  is  an  assembly  of  one  unit  each  of  P?  and  P3. 

Product  2  (P2)  is  composed  of  2  units  of  Froduct  3  (P3) . 

There  are  two  types  of  operation?,  type  1  (T-^) ,  and  type  2  (To). 

One  hour  of  T-^  is  required  for  one  unit  of  P, . 

One  hour  of  T^  plus  one  hour  of  T2  is  required  for  one  unit  of  ?2* 

Two  hours  of  To  are  required  for  one  unit  of  P0. 

The  above  assumptions  are  summarized  in  the  following  table; 


To 

make 

To  make 

To  make 

one 

-  unit  P-^ 

one  unit  P2 

one  unit  Po 

P1  (unit) 

P2  (unit) 

1 

P3  (unit) 

1 

2 

T-,  time  (hours) 

1 

1 

T2  time  (hours) 

1 

2 

24..  Jackson,  J. P..,  "Mathematical  Models  with  Examples  from  Linear 

Pro  jamming.  *'  Industrial  Logistics  Research  F reject  Report  No.  W, 
Los  Angeles,  Univ.  of  California. 
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Further  assumptions: 

3000  hours  per  week  of  T,  are  available. 
4.000  hours  per  week  of  Tp  are  available. 
3000  units  of  P3  are  in  storage  at  the  commencement  of  the 
production  period. 

Net  profit  from  the  sale  of: 

?1  =  m 

P2  s  *  5 

p3  =  e  1 

The  net  increase  in  profit  from  the  sale  of: 
p   -  14  -  5  -  1  -  $8   and  of 

P2  =   5  -  2  (1)   -  13 

Items  produced  during  any  given  week  cannot  enter  into  further 
production  until  the  next  week.  ie.  units  of  P2  made  during  the  first 
week  cannot  be  assembled  into  P,  until  the  second  week. 

Recuirement:  to  operate  for  three  weeks  so  as  to  maximize  profits 

Let: 
x*  a  production  of  P-.  in  week  t:   t  =  1,  2,  3. 
yt     ■      "     »  ?2   in  week  t:   t  =  1,  2,  3. 

zt  -      "      p3  in  week  t;   t  =  1,  2,  3. 

vr1     =     number  of  units  of  P^  stored  during  week  tj       t     =     1,    2,   3. 

vt      =  n  «  11  n      p2  n  11  it       If.  t      s      1,    2,    3. 

wt      _  „  11        n        «     p  n  n  «      ";         t      =     1,    2,    3. 

a*     =    number  of  hours  of  idle  time  of  Tj_  in  week  t;     t     »     1,   2,  3. 
b*     ■     number  of  hours  of  idle   time  of  T2  in  week  t;     t     ■     1,   2,  3. 
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From  the  table  it  can  be  seen  that: 

0  =  number  of  units  of  P^  consumed  in  week  t;  t  »  1,  ?,  3. 

•yt    B        If      !l       it     n    p^ 

*2  "  »  «  »j  t  =  1,  2,  3. 
xt+2yt  =  »  "  "  »  p  «  n  "  ".  t  =  1,  2,  3. 
In  addition: 

0  ■  number  of  units  of  P-^  available  at  beginning  of  week  1. 

0      =  "  n  "  "      p  "  n  n  nun 

2 
3000      ■  n         h  n        tt     p  it  it         n  n         it         it 

u*  +  x*  =  P-,  available  at  beginning  of  week  (t  *  1)  j  t  =  1,  2,  3. 

vt  +  yt  =  p^  ti  n  it  n  "  (t  +  1);  t  =  1,  2,  3. 
v^z*  =  P3  "  "  "  »  "  (t  *  1);  t  =  1,  2,  3. 
From  the  above  restrictions  the  following  equations  can  be  obtained: 

u1  =  0 

v1  +  x1  =  0 

w1  ♦  x1  4-  2V1  =  3000 
I.    u1  =  u1-1  *  x1"1; 

v^  4-  x   «■  v    *■  y^1 

w*  t  x*  «•  2y*  =  w1"1  t  z1"1; 

Availability  of  machine  time  leads  to  the  following  equations: 

II.    xt  4  y1  *  at  =  3000  t  -  1,  2,  3. 

yt  +  2Z1  +  b*  a  4000  t  .  1,  2,  3. 

In  addition: 

III.  x1,  y*,  z*,  u1,  v1,  w*,  a1,  bl,  >  0;      t  =  1,  2,  3. 

Finally  the  net  profit: 

IV.  P  =  S^x1  +  x2  +  x3)  ♦  3(y1  *  y2  *  y3)  *  (zl  *  z2  +  z3) 

Thus,  to  operate  for  three  weeks  so  as  to  maximize  profits  calls 

for  maximizing  IV  subject  to  the  restrictions  imposed  by  I,  II,  and  III. 
The  preceding  problem  by  J.R.Jackson  can  now  be  solved  in  the 


t 

s 

2,   3. 

t 

■ 

2,  3. 

t 

- 

2,3. 

A? 


following  way: 

Substituting  the  values  for  (t)  the  following  equations  are 
obtained : 

u1  »  0 

v1     4  x1  ■  0 


w1  4-  xl 


uJ 


u3  ■  u2 


v*- 

v3 

w2 

v3 

vl 


J1 

v2 


x 
x3 

vl 


2yl  -  3000 

xl 

x2 

v1  f  yl 


*  y2 


2y2 
2y3 


w1  4-  z1 


4  z' 


a1  a  3000 
a2  a  3000 


a3  »  3000 
2a1  ♦  b1  a  ^00 
2z2  +  b2  a  4000 
2Z3  f  b3  ■  ^000 
Since  vr   ■  0  it  may  be  eliminated,  and  since  v*  =  -x^  and  u2  =  x^, 
v  and  u2  may  be  replaced  by  -or*  and  x-*-  respectively.  The  result  is  a 
set  of  12  equations  as  follows: 

1.  w1  4  x1  4-  2^  »  3000 

2.  u3  -  x1  -  x2  =  0 

3.  v2  4-  x2  f  x1  -  y1  ■  0 


A.  v3  +  x3  - 


-  y 


2  =  0 


5.  w2  +  x2  4-  2y2  -  w1  - 


6.  w3  i  x3  4.  2y3  -  w2  - 


,2     - 


a  0 


a 


7.  x1  +  y1  +  a1  -  30OO 

8.  x2  ♦  y2  *  a2  =  3000 

9.  x3  4.  y3  f  a3  =  3000 

10.  y1  4.  2551  •»•  bl  -  4000 

11.  y2  ♦  2*2  f  b2  =  4000 

12.  y3  ♦  2z3  ♦  b3  =  4000 

These  12  equations  are  displayed  in  the  familiar  arrangement  of 
Table  8.  It  appears  that  it  will  be  possible  to  rearrange  the  table 
so  as  to  provide  twelve  linearly  independent  vectors  to  form  the  initial 
basis.  The  below  listed  steps  will  yield  an  equivalent  matrix  which 
will  have  the  desired  form. 

1.  Replace  row  6  by  the  sum  of  row  6  and  row  5. 

2.  Interchange  columns  v3  and  w^. 

3.  Replfcce  row  4  by  the  sum  of  row  4  and  row  3. 

4.  Replace  row  5  by  the  sum  of  row  5  and  row  1. 

5.  Replace  row  6  by  the  sum  of  row  6  and  row  1. 

6.  Interchange  columns  v2  and  w"^. 

7.  Interchange  columns  u3  and  \r. 

Moving  the  basis  to  the  left  side  of  the  table,  adding  the 
functional  coefficients,  the  stub  column,  and  the  F  and  C  -  F  rows 
completes  the  first  set  of  calculations  as  shown  in  Table  9.  Suc- 
ceeding tables  can  then  be  utilized  to  carry  out  the  steps  of  the 
simplex  method. 


Table  8 


Equa- 
tion 

'  Po 

,1 

x2 

*3 

y1 

72 

r5 

i1 

Z2 

Z3 

«3 

T2 

J 

1 

w 

w2 

v^ 

al 

a2 

83 

b1 

b2 

b3 

(1) 

3000 

1 

2 

1 

(2) 

0 

-1 

-1 

1 

(3) 

0 

1 

1 

-1 

1 

U) 

0 

1 

-1 

-1 

1 

a) 

0 

1 

2 

-1 

-1 

1 

(6) 

0 

1 

2 

-1 

• 

-1 

1 

(7) 

3000 

1 

1 

1 

(8) 

3000 

1 

J. 

1 

(9) 

3000 

1 

1 

1 

(10) 

4000 

1 

2 

1 

(11) 

4000 

1 

2 

1 

(1?) 

4000 

1 

2 

1 

Table  9 


c 

8 

8 

8 

3 

3 

3 

1 

1 

1 

Vec 
tor 

P0 

v1 

U3 

V2 

v3 

u2 

w3 

a1 

a2 

a3 

b1 

b2 

b3 

X1 

X2 

x3 

y1 

y2 

T3 

al 

*2 

*3 

V* 

3000 

1 

1 

2 

u3 

1 

-1 

-1 

v2 

1 

1 

1 

-1 

t3 

1 

1 

1 

1 

-1 

-l 

V2 

3000 

1 

1 

1 

2 

2 

-1 

w3 

3000 

1 

1 

1 

1 

2 

2 

2 

-1 

-1 

a1 

3000 

1 

1 

1 

a 

3000 

1 

1 

1. 

a3 

3000 

1 

1 

1 

bl 

4000 

1 

1 

2 

b2 

4000 

1 

1 

2 

b3 

4000 

1 

1 

2 

F 

C  -  F 

8 

8 

8 

3 

3 

3 

1 

1 

1 
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To  show  the  amount  of  computation  reauired  to  solv6  this  problem 
by  an  automatic  computer  recall  that  advancing  from  one  table  to  the 
next  is  in  effect  the  movement  from  one  extreme  point  to  another.  To 
go  from  one  extreme  point  to  another  by  the  following  steps  requires 
approximately  the  number  of  computer  operations  as  indicated  for  each 
step. 

1.  the  calculation  of  the  unit  functional  weights  (C  -  F) . 
286  operations. 

2.  the  inspection  and  comparison  of  the  unit  functional  weights 
and  the  selection  of  the  one  with  the  largest  positive  value. 
This  determines  which  vector  is  to  be  placed  in  the  basis 

36  operations. 

3.  The  selection  of  the  basis  vector  which  is  to  be  replaced  by 
the  incoming  vector.    60  operations. 

I.   the  calculation  of  the  new  value  of  the  functional.  2  opera- 
tions 

5.  The  determination  of  the  column  vectors  in  terms  of  the  new 
basis.    24.2  operations. 

The  total  number  of  operations  required  then  to  advance  from  one 
extreme  point  to  the  next  is  626. 

Computational  experience  indicates  that,  in  general,  the  number 
of  extreme  points  which  must  be  examined  is  approximately  twice  the 
number  of  equations.  In  this  problem  then,  the  total  number  of  computer 
operations  is  (2)  (12)  (626)  a  15024. 

A  Card  Programmed  Calculator  (C.P.C.)  which  is  a  general  purpose 
floating  point  digital  computer  that  can  perform  100  operations  per 
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minute  would  then  take  approximately  two  and  one-half  hours  of  actual 
computer  operation  to  solve   this  problem.     This  computer,   however,   can- 
not store  enough  information  to  allow  continuous  operation.     The  manual 
bookkeeping  required  for  the  removal  and  re-insertion  of  information 
doubles  or  triples  the  time  reouired,  and  this  problem,   therefore, 
would  require  between  five  and  seven  and  one  half  hours  for  solution. 
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CONCLUSION 

Linear  programming  is  thus  a  mathematical  technique  to  find  the 
values  of  non-negative  variables  or  unknowns  which  are  subject  to  linear 
restrictions  and  such  that  the  values  will  optimize  a  linear  functional. 
It  has  been  shown  that  the  linear  programming  problem  can  be  visualized 
as  a  geometric  representation  of  vectors  or  points  in  imaginary  spaces 
of  many  dimensions.  The  solution  of  a  simple  numerical  example  has 
been  demonstrated  using  the  simplex  method. 

The  requirements  that  must  be  met  before  this  technique  can  be 
used  for  practical  problems  have  been  discussed.  It  can  thus  be  seen 
that  linear  programming  has  s  restricted  usefulness  for  Industrial 
Engineering  applications.  As  Joseph  0.  Harrison  says:  * 

The  practical  difficulties  encountered  in  applying  linear  programming 
are  threefold:   (l)  the  expression  of  realistic  objectives  and  con- 
straints in  measurable  terms,  (2)  the  determination  of  suitable  numer- 
ical values  for  coefficients,  and  (3)  the  computational  labor  required 
to  numerically  execute  large  linear  programming  problems. 

Like  any  mathematical  technique,  linear  programming  must  not  be 
considered  a  panacea  for  all  the  problems  of  Industrial  Engineering.  It 
is  restricted  in  the  scope  of  its  application  both  by  the  nature  of  the 
problem  and  by  the  economic  considerations  of  computation. 

The  history  of  the  application  of  linear  programming  to  practical 
problems  covers  too  short  a  period  of  time  and  the  application  has  been 
too  limited  to  permit  a  complete  evaluation  of  the  usefulness  of  this 
technique  to  Industrial  Engineering.  It  can  be  expected,  however,  that 


25.  Harrison,  J.O.,  Jr.,  Linear  Programming  and  Operations  Research. 
Informal  Seminar  in  Operations  Research,  Paper  No.  2;  Baltimore, 
Johns  Hopkins  Univ.,  October  13,  1953,  p.  18. 
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the  considerable  amount  of  study  now  being  devoted  to  both  the  theory 
and  the  practical  aspects  of  linear  programming  will  result  in  the 
overcoming  of  many  of  the  difficulties  limiting  its  usefulness,  and 
that  Industrial  Engineering  will  make  more  and  more  use  of  this  tech- 
nique in  the  future. 
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